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Abstract 

We demonstrate existence of waves localized at the interface of two 
nonlinear periodic media with different coefficients of the cubic non- 
linearity via the one-dimensional Gross-Pitaevsky equation. We call 
these waves the surface gap solitons (SGS). In the case of smooth sym- 
metric periodic potentials, we study analytically bifurcations of SGS's 
from standard gap solitons and determine numerically the maximal 
jump of the nonlinearity coefficient allowing for the SGS existence. We 
show that the maximal jump vanishes near the thresholds of bifurca- 
tions of gap solitons. In the case of continuous potentials with a jump 
in the first derivative at the interface, we develop a homotopy method 
of continuation of SGS families from the solution obtained via gluing of 
parts of the standard gap solitons and study existence of SGS's in the 
photonic band gaps. We explain the termination of the SGS families 
in the interior points of the band gaps from the bifurcation of linear 
bound states in the continuous non-smooth potentials. 



1 Introduction 

We are concerned with localized waves at the interface of two periodic non- 
linear media called surface gap solitons (SGS). One of the first publications 
on optical solitons propagating along material interfaces is jl7] . where the 
interface of a linear and a focusing Kerr nonlinear medium is studied. In 
the last two years relevant publications in the context of nonlinear optics 
have dealt, for instance, with discrete surface solitons in nonlinear waveguide 
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arrays [UllOlfTB]. surface gap solitons at the interface of a uniform and a pe- 
riodic medium with the defocusing cubic nonUnearity [7] and surface vortex 
sohtons at the interface of two periodic media with different mean values of 
the refractive index and with the saturable nonlinearity |5i [6] . One of the 
typical models employed in the theory of gap solitons is the one-dimensional 
nonlinear Schrodinger equation (NLS) with cubic nonlinearity and periodic 
potential called the Gross-Pitaevsky equation. 

We investigate here the existence of surface waves at the interface of two 
media with identical periodic linear parts of the refractive index and with 
different cubic nonlinearities. It is known that for most photonic materials 
a variation in the nonlinear part of the refractive index n2 is necessarily 
accompanied by a larger change in the linear part ng. Nevertheless, certain 
materials exhibit large variations in n2 accompanied by small variations in 
no, see [U [15]. Localized states have been studied theoretically in media 
with constant no and spatially periodic n2 in [3j. 

Each of the two periodic nonlinear media supports at least two families 
of standard gap solitons in every bounded non-empty frequency gap. One 
family is always unstable, while the other can be stable depending on the 
locations of spectral bands and bifurcations of eigenvalues from the band 
edges jl3] . The potentially stable family looks like a single- humped envelope 
soliton with exponential decay and oscillations near the central peak. Multi- 
humped envelope solitons may also exist in such periodic nonlinear media 
but we shall focus herein on existence of a single-humped solution localized 
near the interface between the two periodic nonlinear media. 

The paper is organized as follows. Section [2] reviews Floquet theory 
for the governing Gross-Pitaevsky equation and summarizes the results on 
existence of gap solitons. In Section [3] we study the existence of SGS's for a 
smooth symmetric periodic potential function and find the maximal allowed 
jump in the nonlinearity coefficient between the two media for existence 
of SGS's. Section [5] discusses bifurcations and existence of SGS's for a 
continuous potential function with a derivative jump at the nonlinearity 
interface. Section [5] concludes the paper with conjectures on stability of 
SGS's. 
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2 Background: Floquet theory and gap soliton ex- 
istence 

We consider the one-dimensional periodic cubic Schrodinger equation in the 
form 

iut = -Uxx + V{x)u-r{x)\u\'^u, xeR, t>0, (2.1) 

where x and t are the spatial and temporal variable respectively, V{x) is a 
real, continuous and d-periodic potential, and T{x) = T± for ztx > is a 
real nonlinearity coefficient with constants r_|_ and r_. The positive values 
of T(x) corresponds to the focusing nonlinearity and the negative values of 
r(x) to defocusing nonlinearity. 

We are interested in the existence of stationary solutions of (j2.1|) local- 
ized near the interface at x = and having the form 

t) = e"*'^V(2;) s.t. 0:M^M, (/> ^ as |x| ^ oo. (2.2) 

The function (j){x) has to satisfy the second-order non-autonomous ODE 

-(p" -uj(l) + V{x)(p-r{x)(l)^ = 0, (2.3) 

which can be cast in the Hamiltonian form with the Hamiltonian function 

m = \ w? + ^^0' - v{x)c^''] + \nx)c^\ (2.4) 

Since r(x) is discontinuous at x = 0, (t){x) is a weak solution of the ODE 
(j2.3p in G C^(M+ U M_), such that the second derivative 4>"{x) may have 
a jump at X = 0. The continuously differentiable solution cj) G C^(M) is a 
critical point of the energy functional 

E.[^] = \ [ [-\<P'\' + u;\^\'-VixM'']dx + ^ [ rixM\Ux, 

such that the first variation of E'^[(f)] recovers the ODE (j2.3p . 

Replacing t by z in (j2.1|) and (j2.2p . the x— localized solution n(x,z) can 
be viewed as a spatial soliton propagating along the direction z and localized 
in the transverse direction x (localization in the third spatial direction y is 
assumed to be achieved via total internal reflection). The parameter uj plays 
a role of the propagation constant. As we show below, the localized solutions 
of the ODE (j2.3p decay exponentially as |x| — > oo only if a; belongs to the 
frequency gaps in the continuous spectra of the operator L := —dxx + V{x) 
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called the photonic band gaps. To do so, we recall the basic Floquet theory 
(see [21 E]) for the Hill's equation 

L^l,{x) = -ip'^x) + V{x)'ip{x) = ujij{x), x e M. (2.5) 

The bounded solutions 'ip{x) of the Hill's equation (j2.5p are usually called 
Bloch functions. Given a real, continuous and (i-periodic potential V{x), 
bounded solutions tplx) exist for w in a union of (possibly disjoint) spectral 
bands from 

S := [ujQjUJi] U [W2,t^3] U [uJi,uj5] U . . . , 

where W2n-2 < ^2n-i < ^2n, n € N and — > cx) as n — > oo. The set 
S represents the complete (purely continuous) spectrum of the operator L 
[2]. We shall assume for simplicity that all spectral bands are disjoint with 
'^2n-i < ^2n, n gN, such that all finite frequency gaps are non-empty. 

For a fixed lo in the interior point of S, both fundamental solutions of the 
second-order ODE (j2.5p are quasi-periodic in x and have the representation 
ip = p±{x)e^^''^ , where p±{x) = p±{x + d) and k G [O, ^] . The parameter 
k parameterizes the frequency parameter lo, such that we shall use notation 
= ^2n,2n+i{k) for the spectral band in w E [u2n,i^2n+i]- If the n-th band 
is separated from the (n + l)-th band (i.e. ui2n-i < ^2n and W2n+i < ^2n+2), 
then uj2n 2n+i(^) — ^ cJ2n 2n+i(^) 7^ at the end points A; = and k = ^ 

When u = uJn, one of the solutions ip = ipnix) is either d-periodic (cor- 
responding to A; = 0) or d-antiperiodic (corresponding to A; = ^) and the 
other fundamental solution ^l^{x) grows linearly in x. For a fixed a; € M \ S 
the two fundamental solutions of ()2.5p grow exponentially either in x or —x 
and have the representation = u±{x)e^'^^ , where u±{x) is either periodic 
or anti-periodic and k = k{uj) € M_|_. The functions u±{x) are periodic 
(anti-periodic if the bounded solutions Tpn{x) are periodic (anti-periodic) at 
the band edges W2n-i and u!2n, which surround the band gap. 

Suppose that (j){x) is a localized solution of the ODE (12. 3|) . It is then 
obvious from the linearized analysis that the solution (p{x) decays exponen- 
tially as |x| oo only if w G M \ S. It was shown under fairly general 
assumptions (see [T^ and references therein) that the families of gap soli- 
tons of the ODE (j2.3p with constant coefficient T{x) = Tq undertake a local 
bifurcation from all points oj = a'2m; m > to the left if Fq > and from all 
points LJ = uj2m+i, m > to the right if Fq < (the term local bifurcation 
means that — > as c<j ^ w„). This conjecture was rigorously proved 
in where existence of exponentially decaying gap solitons in H^(M.) was 
confirmed in every finite frequency gap lv € {L02m-i,^2m), m S N and in 
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the semi-infinite frequency gap to < ujq lor Tq > 0. We use this result but 
simplify our consideration by working with the class of symmetric potentials 
V = Vo{x), where Vq{—x) = Vo{x) on x G M. In particular, we shall perform 
numerical computations with 



which has a minimum at x = 0, i.e. at our interface location. The spectral 
bands and gaps of Vo{x) are approximated numerically from the Hill's equa- 
tion (j2.5p . For instance, the first five band edges of the potential ()2.6p are 
located as follows 

ojQ « 0.283, wi 0.291, 002 ~ 0.747, ^3 ^ 0.843, 0^4 ~ 1.057. 

As partly seen in Fig. 1 of ^13j, for Vq{x) = sin^(7rx/10) the Bloch functions 
■0 = ipnix) at the first eight band edges to = tOn, n E {0, 1, . . . , 7} have the 
following symmetry properties: 



Clearly, all ipnix) must be even or odd since the Hill's equation (j2.5p is 
symmetric with respect to reflection x 1— > —x and admits only one linearly 
independent bounded eigenfunction = ipnix) at a; = iVn- Nevertheless, 
for other even potentials Vo(x) the ordering between even and odd Bloch 
functions can be different than in (j2.7p and the only statement about this 
ordering that is valid for general even Vo(x) can be deduced from Theorem 
3.1.1 in [2j, which says that two subsequent even eigenfunctions have to be 
followed by an odd one and vice versa. We refrain from such a discussion 
here concentrating on the potential ()2.6p . for which symmetries (|2.7p hold. 

In contrast to the ordering between even and odd Bloch functions, the 
ordering between d— periodic and d— antiperiodic Bloch functions is unique 
for all d— periodic Vo{x) (not necessarily even). By the trace of the mon- 
odromy matrix [2j, the periodic functions ipn{x) correspond to the set n G 5"+ 
with S+ = {0,3,4,7,8,...} = {0,4A; - l,4/c;A: G N} and the anti-periodic 
functions Tpn{x) correspond to the set n G S- with 5_ = {1,2,5,6,...} = 
{4A;-3,4A:-2;feGN}. 

By Sturm's Theorem (see Theorem 3.1.2 in [2]) for general d— periodic 
Vb(x) the periodic functions tpnix) with n G 5+ have no nodes on x G [0, d) 
for n = 0, two nodes for n = 3, 4, four nodes for n = 7, 8, etc. Similarly, the 




(2.6) 



0„(-x) = ipn{x), n G {0, 1,4,5}, 
ipni-x) = -Tpnix), u G {2,3,6,7}. 



(2.7) 
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anti-periodic functions ipn{x) with n ^ S- have one node on x G [0, d) for 
n = 1,2, three nodes for n = 5, 6, etc. 

Altogether, this set of facts about Bloch functions at the first eight band 
edges for the potential (j2.6p is summarized in Tabled) 
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-1 



Table 1: Properties of the Bloch functions ipn{x) and gap soliton bifurcations 
at the first eight band edges of the even potential Vq{x) = sin^(7rx/10). 

Let (t>o{x) be a single-humped solution of the ODE (|2.3p with T{x) = Tq 
and V{x) = Vo{x) = sin^(7rx/10) which bifurcates from the band edge u! = 
u!n- By the local bifurcation theory [13], it inherits the symmetry properties 
()2.7|) of the Bloch function Tpnix). Therefore, (j)o{—x) = (poix) for branches 
of gap solitons to the left of w„ with n = {0, 4} (for Fq > 0) and to the right 
of LOn with n = {1,5} (for Tq < 0), while (poi—x) = — 0o(x) for branches of 
gap solitons to the left of w„ with n = {2, 6} (for Tq > 0) and to the right of 
LOn with n = {3,7} (for Tq < 0). See Figs. 2-3 in [13j for gap solitons (j)o{x) 
in the potential (|2.6p . 

In this paper, we shall consider existence of surface gap solitons in the 
ODE (j2.3p with piecewise constant coefhcient r(x) = T± for ±x > and 
potential V{x) of the following two classes: 

{i)V = Vo{x), {u)V = Vo{x-6)xi-oo,o) + Vo{x + 6)x[o,oo), (2.8) 

where X[a,b] = 1 on x G [a, b] and zero otherwise, while < 5 < d. Here 
Vb(x) is a smooth, even, d-periodic function on x € M. We note that V{x) 
in (ii) is continuous and even on x € M but smooth and periodic only on 
each ibx > 0. 

One can develop a general shooting method for numerical approxima- 
tions of SGS's from the condition that a localized solution (p{x) of the 
second-order ODE (12. 3p with to E {uj2m-i,^2m), m G Z decays to zero 
at infinity according to two fundamental solutions p±(x)e^'^^ as x — > ±oo, 
where k = is a positive number. Solving the ODE (j2.3p with r(x) = 
for a general initial value <^(0) and (j)'{0) to x > and the same ODE with 
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r(x) = r_ to x < 0, one can construct a continuously differentiable solution 
(j){x) on x G M which decays to zero as x — > ±00 if and only if the projec- 
tions to the growing fundamental solutions p±{x)e^'^^ are zero at infinity. 
The system of two constraints for two initial values constitute a well-posed 
problem of numerical analysis. However, practical implementations of this 
algorithm are unclear as the shooting method may depend sensitively on 
starting approximations of the initial value and may require long computa- 
tions time to search through all appropriate initial values. In addition, the 
ODE solvers of the shooting method may develop numerical instabilities in 
approximations of growing solutions. 

Due to these reasons, we shall develop an alternative view on numerical 
approximations of SGS's, starting with local bifurcation analysis and using 
the homotopy continuation method to trace the solution families along pa- 
rameters uj, T± and 5. Using these analytical and numerical results, we have 
obtained the following main results of the article. 

1. We prove analytically that any gap soliton for r+ = r_ can be contin- 
ued to the SGS for sufficiently small IF^ — r_| under a non-degeneracy 
assumption. 

2. We prove analytically that the maximal difference |r_|_ — r_| leading to 
SGS existence converges to when to approaches the band edge which 
features the local bifurcation of a gap soliton. 

3. Surface gap solitons are computed numerically when the potential 
V{x) is given by (I2.8p (i). and the maximal |r+ — r_| allowing their 
existence is found. Our numerical results confirm the analytical results 
(1-2) above. 

4. Existence of SGS's for V{x) in (ESKfi) with r+ > and r_ < 
is studied. We show numerically that local bifurcations may occur 
from a countable set of points in the parameter domain {uj,6) G 
(w2m-l,W2m) X (0,d), m G N. 

5. We compute numerically the points of local bifurcation of SGS's for 
the potential (12. Sh (ii) and use the homotopy continuation of the bifur- 
cating solution. As a result, we show that the family of SGS's exists 
typically in a subset of the plane (w, 5). 

6. We show analytically that the termination of families of SGS's for the 
potential (j2.8p (ii) is related to existence of linear bound states for the 
non-smooth potential. 
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Results (1-3) are reported in Section [3] and results (4-6) are described in 
Section [H 

3 Bifurcations of surface gap solitons for smooth 
potentials 

In this section we study continuation of SGS's from gap solitons existing for 
F-i- = r_ in the case of a smooth potential function V{x). A prototypical 
example of such potential is the symmetric function Vo(x) in (|2.8|) (i). 

3.1 Existence of Bifurcations from Gap Solitons 

Let 7 = (r+ + r_)/2 and = (r+ - r_)/2. Then, the ODE ([23]) can be 
rewritten in the form 

F{(1), u) = -(j)" -uj(l3 + y(x)0 --i(t)^ -V sign(x)(/)^ = 0, (3.1) 

where -P((/), z^) : i?^(M) x M i-^ //^^(M) is a nonlinear operator acting on a 
function ^(x) in space (j) £ //"^(M) and parameter G M. 

We assume that there exists a solution 4)q{x) G H^(U.) for o; G M \ S and 
some 7 and V{x), such that F(0o, 0) = 0. The Jacobian D^pF {(f>o,0) is given 
by the Schrodinger operator C : //^(M) i-^^ L^(R), where 

C = -dl-uj + V{x) -3'y(f)l{x). (3.2) 

Since w G M \ S, we have (j)Q{x) — > exponentially fast as |x| —i- oo, such 
that the term —3^(j)Q{x) is a relatively compact perturbation to the un- 
bounded operator L — uj, where L = —9^ -|- V{x). By a standard argument 
(see Corollary 2 in Section XIII. 4 in [H]), the essential spectrum of C and 
(L — uj) coincide. Since uj G M \ S, the zero point is isolated from the es- 
sential spectrum of >C. If we further assume that C has the trivial kernel in 
Fi(R), then C is invertible on L^{R). Since the translational invariance is 
broken if V{x) ^ 0, C generally has the trivial kernel, unless a bifurcation 
of branches of gap solitons occur. By the standard analysis based on the 
Implicit Function Theorem, there exists a unique smooth continuation of 
4>uix) from (/i>o(x) in //^(M) for sufficiently small such that F{(f)y,v) = 
and (t)v{x) 4'o{x) in H^{M.) as u 0. 

In other words, we have proved above that if a gap soliton exists for 
r_|_ = r_ and a; G M \ E and the linearized operator C is non-degenerate, 
then the gap soliton is uniquely continued into the SGS for small non-zero 
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|r+ — r_|. We confirm this prediction via numerical analysis of the ODE 
(|2.3p with V{x) in (|2.8p (i) for lo taken in the semi-infinite band gap and the 
first two finite gaps. Numerical approximations of (poix) for r+ = r_ are 
obtained from the Newton-Raphson iterations and the homotopy continu- 
ation method. The initial guess for the Newton's iteration is taken from 
an asymptotic expansion leading to the NLS approximation [13] when co is 
close to the local bifurcation threshold a;„. After a successful convergence 
for one such co we use a standard homotopy continuation and generate a 
family of gap solitons (poi^) parameterized by uj. The discretization of the 
ODE (j2.3p is based on a fourth order central difference approximation of 
dxx on a truncated domain with zero Dirichlet boundary conditions. 

3.2 Numerical Computations of Surface Gap Solitons 

We now proceed to construct SGS's, i.e. solutions (j){x) of the second-order 
ODE (12. 3p with r+ 7^ r_. When (j)o{x) is obtained for a given value of 
LO, we can apply the numerical homotopy continuation of the solution by 
deviating r_ from r+. At each step, the SGS (j){x) is thus found via Newton's 
iterations. The final value of r_, up to which the iteration converges, is 
denoted by E^,. 

Fig. [T] shows the values of for E_|_ = +1 (a) and E-|- = —1 (b). The 
computational tolerance in E* is 0.006 inside the band gaps and 0.002 near 
the band edges. In the case E_|_ = 1, local bifurcations of small-amplitude 
gap solitons occur from the lower band edges. Eig. [T]^a) shows that the SGS's 
exist in the semi- infinite gap, as well as in the first two frequency gaps. In 
the case E+ = — 1, local bifurcations of gap solitons occur from the upper 
band edges. Eig. [IJb) shows that the SGS's exist in the first and second 
frequency gaps. The two insets of Eig. W[a) show that E* decreases fast 
as UJ moves away from the edge of the first band and that the convergence 
E^, I 1 as LiJ I Wo is smooth. We further see from Eig. [T]that the interval of 
existence shrinks as uj approaches the value uJn for any band edge, where gap 
solitons undertake a local bifurcation. In addition, the interval of existence 
is extremely large in the semi-infinite gap (— oo,a;o), but it becomes narrow 
in the finite gaps {uJ2m-i,^^2m) for m > 1. 

Eor comparison, the family of SGS's in the gap {uji,uj2) exists for —0.24 < 
E,, < 1 in the case E_|_ = +1 and —1 < E^, < 0.47 in the case E_|_ = — 1. 
The family of SGS's in the gap (013, a;4) exists in a very narrow region of 
0.92 < E* < 1 in the case E+ = +1 and in a bigger interval — 1 < E^, < 0.37 
in the case E+ = — 1 (similarly to that in the first gap). 

Fig. [2] shows profiles of SGS's which correspond to the twelve points 
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(a) Upper edge SGS families: r =1 




(b) Lower edge SGS famiiies: r = -1 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 



% 0)^ W (1)3 (H^ 

Figure 1: The values of for SGS's originating from symmetric GS families 
of the first three frequency gaps of V{x) = sin^(7rx/10). In (a) the upper 
inset zooms in and the lower inset zooms out on the graph in the semiinfinite 
gap. The points A — L are referenced in Fig. [2l 

labeled A — L in Fig. [TJ The full lines correspond to the gap solitons from 
which the homotopy in F_ is started (i.e. points A,D,G and J). Clearly, 
the total power and maximum amplitude of the SGS increase as |F_|_ — F_| 
increases. Also notice that the profiles become more concentrated on the 
half X > in the case r+ = +1 [see Fig. [2] (a-b)] and on the half x < in 
the case r+ = —1 [see Fig. [2](c-d)] as |F_|_ — F_| increases. This is in accord 
with the law of refraction: when F_|_ = +1 and F_ decreases from 1, the half 
X > becomes relatively more focusing and therefore attracts more energy 
of the soliton, while when F_|_ = —1 and r_ increases from —1, the situation 
is opposite. 

3.3 Asymptotic Analysis near Gap Soliton Bifurcation Points 

We shall explain now why the existence interval shrinks to zero when uj 
approaches the value ujn where a local bifurcation of gap solitons occurs. As 
UJ — > Un, we have ||(/'o||l°° and C —i- {L—uJn)- Since the operator (L—LOn) 
is not invertible, the Implicit Function Theorem can not be used and the 
solution (poi-c) can not be continued beyond = 0. In order to give a more 
precise explanation of this phenomenon, we adopt the NLS approximation 
for local bifurcation of gap solitons from [13j (see also review in |12j). In 
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Figure 2: The profiles of SGS's corresponding to the points A — L in Fig. [H 
Values of oj are A- C: 0.26, D-F: 0.55, G-I: 0.6, J - L: 0.96. Values of 
r_ are A: 1, B: -3.9, C: -15.3, L»: 1, E: 0.38, F: -0.235, G: -1, -0.45, 
/: 0.002, J: -1, iC: -0.4, L: 0.164. 

particular, we consider an asymptotic solution to the ODE ()2.3p : 

(3 3) 

</.(x) = £^(X)V'n(2;) + A'(X)V;n(x) + X) + 

where X = e x, e « 1, the function A{X) and parameter are defined 
below; ipn and ipn are the d-periodic (or d-antiperiodic) Bloch function and 
generalized Bloch functions respectively of the Hill's equation (|2.5p for lo = 
ujn, such that 

(L - ujn)i>n = 0, (L - a;„)^„ = 2<. (3.4) 

The correction term <j)'^'^\x,X) at C'(e'^) solves the non-homogeneous prob- 
lem 

{L - uJn)(t>^^'' = ^AiJn + A"iJn + 2A" ^'^ + T{X)A^^l. (3.5) 

To ensure boundedness of (l)^^\x,X) with respect to the variable x, and, 
hence, legitimacy of the expansion (j3.3p . one has to apply the Fredholm 
alternative which imposes the orthogonality condition of the right-hand-side 
of (j3.5p with respect to tpnix) on x € [0,d]. The orthogonality condition is 
written as follows 

nA + fiA" + pT{X)A^ = 0, (3.6) 
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where 

(V'n,V'n)' ('0n,'0n)' 

and we have used the standard inner product (•, •) over one period x € 
[0, d]. It is shown in p3] that /x = \uJ2n 2n+i(^) with either A; = or A; = 
at the point w = u;„, where 0J2n,2n+i{k) is the dispersion relation between 
u € [tJ2n,W2n+i] and k G [0, §]. 

Due to the nature of the nonUnearity interface, the function T{X) is the 
same as r(x), i.e. T{X) = T± for itX > 0. We shall prove that no localized 
solution of the ODE (|3.6|) exists under the condition r_ 7^ r_(_. Indeed, 
consider the Hamiltonian of the ODE (j3.6p : 

H[A] = \ {i,{A^f + + \pT{X)A\ (3.7) 

If A{X) solves the ODE ([M]), then 

= \pV\X)A\X) = \p{^T^-VS)b{X)A\X\ 

where b{X) is the Dirac delta-function. If A{X) is a localized solution on 
X G M, then the integration on X € M gives the constraint 

0= hrn H\A{X)\- lim H\A{X)\ = -p{T ^ -V J)A\^), 

since ll{A{X)\ ^ if A{X),N{X) ^ as jX| ^ 00. Therefore, ^(0) = 
if r+ 7^ r„. Consider now i?[yl(X)] on X > 0. It is clear from the decaying 
conditions as X — > 00 that i?[74(X)] = const = 0, which together with the 
fact ^(0) = leads to = \miH{A[X)\ = ^p\A'{0)\'^, such that A'{0) = 0. 

The only solution of the ODE ([3:6]) with A{0) = A'{0) = is the zero 
solution A{X) = 0. 

If r+ = r_ = To and sign(/i) = sign(pro) = -sign(J^), the ODE dM} 
has the standard sech-soliton decaying as \X\ 00. However, the result 
above shows that the sech-soliton with = r_ cannot be homotopically 
continued to a decaying solution of (|3.6p for r+ 7^ r_. This proves that 
r^, — >■ r_|_ as a; — > a;„ where u;„ is a local bifurcation value. 

4 Bifurcations of surface gap solitons for nons- 
mooth potentials 

In this section, we study local bifurcations of solutions of the ODE (j2.3p 
when V{x) is a continuous function with the jump in the first derivative 
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at the nonlinearity interface. The prototypical example of such potentials 
is given by (|2.8p (ii). where Vo(x) is an even potential (in our numerical 
computations we use Vq from (j2.6p ). We shall consider the existence of 
SGS's under the normalization r_|_ = — r_ = +1. 

4.1 Surface Gap Soliton Numerical Construction via Gluing 

The point (5*, w*) in the parameter domain 5 G (0, d) and to G {i02m~i,'-^2m), 
m G N is defined to be a point of a local bifurcation of SGS's according to 
the following 2-step algorithm. 

(i) Construction of Continuous Solutions 

Let (l)±{x;uj) denote the family of single-humped gap solitons parame- 
terized by a; G (ct'2m-i) '^2m) and centered at x = corresponding to the 
equation (j2.3|) with T{x) = F-i- respectively. These families bifurcate from 
the points u = uj2m for r+ > and uj = uj2m-i for r_ < 0. In order to 
find continuous solutions, we now study for each fixed 5 G (0, (i) the two 
functions 

/a(w) = </>_(-5;cj) - (t)+{5;u), fsiuj) = ())^{-6;uj) + (p+{6;uj) 

and find their zeros denoted by uja,b = ^a,b{S), respectively. For each 5 
existence of zeros of either /a(i^) or fsi^) is guaranteed by continuity of 
(j)± as functions of u) and by the fact that <j)-{—5;L02m-i) = 4'+i^'^^2m) = 
and (j)^{—5;u)2m) 7^ 0, 4)+{6]L02m-i) 0- Moreover, several zeros of these 
functions may occur for the same 6. 

When a zero coa{S) or ujb{S) is found, a (5-parameterized family of con- 
tinuous solutions (j)Aix;5) or (j)B{x]5), respectively, is constructed by gluing 
two individual gap solitons 

(l)Aix;5) = (j)-{x - 6;uja)x{-oo,o) + (f^+ix + S;u;a)X[o,oo) 
(t)B{x]S) = </>_(x - (5;c^b)x{-oo,o) - <t'+{x + 5;^b)X[o,oo) 

The functions 4'A,Bix', 5) decay as |x| oo and are smooth in x everywhere 
except at the nonlinearity interface x = 0, where they generally have a jump 
in the first derivative. 

Note that it is important to consider both (pA and (pB due to the sign 
invariance of the ODE (12. Sp . Each sign produces a branch of continuous 
solutions of the ODE (lOl) . 

Figs. [3] (a) andm (a) present the numerically computed cja,_b('^) in the 
gaps {loi,U2) and (0:3,104) respectively. The lack of smoothness in the curves 



13 



in these figures is due to an insufficient resolution in the search algorithm 
and can be corrected with a finer resolution. Note that when a;yi,_B(<5) is 
multiple valued as seen in Fig. [3] (a) and [H (a) , we may have several decay- 
ing solutions (j)Aix) and/or (I)b{x) for the same 5. 

(ii) Construction of Surface Gap Solitons 

Next, we search for continuously differentiable solutions within the above 
family (l)A,Bix; 6). To ensure the continuity of the first derivative of (j){x; 6) 
at X = 0, we search for zeros of the two functions 

g^[5) = - (A+(<5;u;a), gsiS) = (I)'_{-S;u;b) + cI)'+{5;ujb)- 

If a zero of either (7^(5) or gsi^), denoted by 5*, exists, then the function 
(j)A{x;6*) or (^5(2;; 5*), respectively, in (14. ip has a continuous first derivative 
across the point x = 0. Figs. [3] (b) andH] (b) present the numerical results 
on computing 5^,. The labelled intersection points O, P,Q, R, S and T cor- 
respond to zeros of gA,BiS)- They are found as intersection points of solid 
and dashed curves of the same color. The solid black line shows the plot 
of (j)'_^_{6;uJA) and the dashed black line shows (/)'_{— S;uj a)- Similarly, the 
solid gray line plots —(P'_^_{6;u;b) and the dashed gray line plots (/)'_{— 6 -,10 b)- 
Therefore, an intersection of a solid black and a dashed black line (points 
0,Q,S) gives zeros 6^ of gAi^) and, thus, a SGS (pAix;6^). Similarly, 
an intersection of a solid gray and a dashed gray line (points P, R, T) gives 
zeros 5* of gB{^) and, thus, a SGS (/>b(x; (5^,). 

Table [2] shows the approximate computed values of (5* and corresponding 
(J* = uja,b{^*) at the points O — T for branches A,B of solutions given by 
(|4.ip . Note that additional points ((5*,u;*) can be obtained by generalizing 
the above functions fA,B and gA,B to 

fjAi^) = 4'-{-Ud+6);u;)-(l)+{jd+6;u;), fjei^) = (l)-{-{jd+6);uj)+(l)+{jd+d;uj) 
and 

9jAi^) = (p'~i-{jd+Sy,u;A)-(p'+{jd+6;ujA), gjsi^) = (t>'~{-{jd+5);uJB)+(t>'+{id+5;u)B) 

for j E {1,2,...} with V still defined as in (j2.8p fiil. Non-trivial points 
(u;*,(5#) may exist for any such j. For illustration we have found one such 
point for j = 1. The computed value is (w*,(^*) (0.73,7.33) and the 
resulting SGS corresponds to the point Z in Fig. [5^a). Such additional 
solutions are SGSs of smaller amplitude compared to those for j = 0. 
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Figure 3: Two-step search for in the gap {u^x^io-})- (a) Result 

of step (i) - parametrization of the famihes of continuous solutions ()4.ip : 
black line ijJa{S)^ gray line ijJb{S)\ (b) step (ii) - search for b^: solid black 
4>\(b\ijJA)^ dashed black (/)'_(— 5; cja), solid gray — </>'^((5; cj^) and dashed gray 
(j)'_{—5;u>B)- Labeled points correspond to SGS's. 
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point 


O 


P 


Q 


R 


S 


T 


branch of solution 


A 


B 


A 


B 


A 


B 




0.58 


0.70 


0.94 


0.97 


1.03 


1.03 


5* 


1.54 


9.66 


0.78 


3.24 


7.97 


9.57 



Table 2: Bifurcation points for surface gap solitons in the domain uj € 
{uji,uj2) U (1^3, u;4) and 6 G (0, d) for branches of solutions given by (j4.1|) . 



4.2 Numerical Homotopy Continuation of SGSs 

Assuming the existence of a point (lj*, d*), we have constructed the SGS of 
the ODE (j2.3|) . where the potential function is given by (j2.8p (ii) and 

(u, 5) = (u*, 5*). The surface gap soliton denoted as is represented by 

one of the functions in (|4.1|) with {to, 5) = {u)^,6^). Each of these solutions 
can be used as a starting point for a numerical homotopy continuation to 
generate a family of SGS's parameterized by w C (<^2m-ii ^^2™) for a given 
value of 6 = 5^:. Similarly, for a fixed w = a family parameterized hj 5 C 
(0, d) can be constructed. Under the same assumption that the operator C = 
—d"^ — cj* + V{x) — 3r{x)(j)l{x) is invertible, the Implicit Function Theorem 
implies that there exists a unique smooth continuation of the particular 
solution (j)*{x) to the family of solutions along parameters oj and 6. 

We restrict our numerical studies to the continuation in lo. Numerical 
results of such continuation from the SGS's at points O — T are shown in 
Figs. [5] (a) and[6](a). The curves plot the total soliton power as a 

function of frequency uj for fixed 6 = 6^. Note that each curve corresponds 
to a different value of 6^ and hence a different potential V{x). The values 
of (5* can be read in Table [2l Termination of a continuation curve is defined 
when the total power of the soliton becomes zero or when Newton iteration 
convergence fails. As the figures show, the latter case is always accompanied 
by the slope of the continuation curve becoming infinite suggesting a viola- 
tion of the implicit function theorem assumptions. The former termination 
case is studied in the following subsection. 

4.3 Analysis of Termination Points of Surface Gap Solitons 

We shall now consider the termination points of the solution families plotted 
in Figs. [5] (a) and [6] (a) where the soliton power becomes zero. The points 
are labeled Ol — 7l and their corresponding 5 and uj are given in Table [31 
The termination points are expected to be related to existence of non- 
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Figure 5: (a) SGS continuation curves, total power versus frequency, in the 
gap (wi, ci;2)- Labeled points O, P correspond to those in Fig. [3] (b). Point Z 
is discussed in Sec. 14.11 Points Ol,Pl and are SGS termination points, 
(b) Point spectrum of the linear Schrodinger operator inside {loi,uj2) for all 
6 G (0, d). Full/dashed lines: eigenvalues with even/odd eigenfunctions. 




(0 Ql T|S|^ (0 

Figure 6: (a) SGS continuation curves, total power versus frequency, in 
the gap {u)3,uJi). Labeled points Q — T correspond to those in Fig. d] (b). 
Points Ql — Tl are SGS termination points, (b) Point spectrum of the 
linear Schrodinger operator inside {ijj^,uj4) for all 5 E (0, d). Full/dashed 
lines: eigenvalues with even/odd eigenfunctions. 
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point 


Ol 


Pl 


Ql 


Rl 


Sl 


Tl 


6 


1.54 


9.66 


0.78 


3.24 


7.97 


9.57 




0.65 


0.70 


0.98 


1.02 


1.03 


1.03 



Table 3: Termination points for the six SGS families in Figs. [5] (a) and [6] 
(a). 



trivial bound states in the (point) spectrum of the Schrodinger operator for 
the same potential V{x)^ i.e. with exponentially decaying solutions of the 
linear ODE 

- V" -0Ji^ + V{x)i: = 0, V : IK ^ IR, (4.2) 

for V{x) in ()2.8p (ii) and to G The point spectrum is nonempty due to 

the singularity of V{x) at j; = 0. 

4.3.1 Numerical Results 

Results of numerical computations of the point spectrum contained in the 
first two finite gaps {uji,uj2) and (tt'3,ci;4) are shown in Figs. [5](b) and[6](b) 
for all values 5 £ [0,d]. The eigenfunctions ip are either even (full lines) or 
odd (dashed lines). For the six values of 6 corresponding to the SGS families 
in Figs. [5] (a) and [6] (a) the eigenvalues are marked by black dots and are in 
perfect agreement with the values of oj at the termination points Ol — Tl- 
The symmetry (even/odd) of the bound states at Ol — Tl also matches 
that of the eigenfunctions at the marked points in the point spectrum. The 
eigenvalue curve originating as well as ending at uj2 in Fig. [5] (b) corresponds 
to the termination point Zl of the SGS family for j = 1 in Fig. [5] (a). The 
termination point Zl for the same value of 5 is shown by a triangle. 

4.3.2 Bifurcation Analysis for \5\ Small 

In the remaining part we consider bifurcations of point spectrum of the 
Schrodinger operator from the band edges for small values of (or, due to 
the d periodicity of V, equivalently for 6 near from above and near d from 
below). This analysis will prove the existence of the spectral curves near 
6 = and 5 = 10 in Figs. [5] (b) and [6] (b), i.e. the existence of curves with 
points Ol and Ql locally to (5 = and the curves with points Pl and Tl 
locally to 6 = 10. 
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In order to construct solutions of the spectral problem (j4.2p . we first 
consider exponentially decaying solutions of the ODE on the half-line 

-V'+-wV+ + Vb(x + (5)'(/'+ =0, 

By using the fundamental solution of the Hill's equation (j2.5p . we can express 
ip+{x) in the form ijj^ = e~'^^U-{x + 5), where U-{x) are periodic or anti- 
periodic bounded solutions of the Hill's equation (|2.5p with V{x) = Vq{x), 
i.e. solutions of (j2.5p at a band edge oj = ojn- 

As V{x) is even, the function ^+{x) admits a symmetric (even) reflection 
about X = if ^^+(0) = 0, which is equivalent to the condition 

Gi((5,k) = n'„((5) - KU-{5) = 0, 

and it admits an anti-symmetric (odd) reflection about x = if ^+(0) = 0, 
which is equivalent to the condition 

G2(5,k) = n_(5) = 0. 

Since eigenvalues of the spectral problem (|4.2p are simple and the eigenfunc- 
tions are either even or odd, all eigenvalues of the spectral problem (j4.2p in 
the band gaps oj G R\S are defined by zeros of the functions Gi{5,k) and 
G2{S,k) in K for a given value of 5, where k > and the values of k are 
related to the values of to in the band gaps. Both functions Gi^2 are ana- 
lytic in 5 G M and periodic with period d. Both functions admit analytic 
continuation in the parameter k € ]R+ [8j. 

If (5 = 0, the only zeros of Gi{5,k) and G2{6,k) occur at k = 0, i.e. 
at the band edges to = ujn- Indeed, if Gi(0,k) = 0, then ip'^{0) = 0, such 
that ip+{x) = ipn{x) is an even function on x € M. However, ip+{x) decays 
exponentially as x — > oo and grows exponentially as x ^ — oo if k > 0. 
Therefore, Gi(0,k) = is equivalent to k = 0. A similar argument works 
for G2(0,k) = 0. 

(i) Bifurcation of Even Eigenf unctions 

Let us first consider the zeros of Gi{5, k). Computing the derivatives of 
Gi{5,k) in 5 and k at {5,k) = (0,0), we obtain 

d5GiiO,0) = u'L{0) = i^';,{0) = {VoiO)-u;n)MO), 

a,Gi(o,o) = -<(o)-Vn(o), 

where ipn is the generalized Bloch function, see (|3.4p . The fact = 
— is clear from differentiation of (j2.5p with respect to k. 
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It is found in [12] that 

D{x) = ijn{x)i)'^{x) - V'^i(a;)V'n(2;) + ijl{x) 
is constant in x, i.e. D{x) = D{0), and that 

^(0) = ^<-l,2nW(V'n,Vn), (4-3) 

where either /c = or /c = ^ at the bifurcation point cj = Since ip'ni^) — 
0, -D(O) = '^pn{0){^p'n{0) + ipniO)) and the leading-order approximation for 
the root of Gi{6, k) near {6, k) = (0, 0) is given by 



</'„(0)(Vo(0)-i.v) ' *5(0)(l'o(0)-^„) 

where V'n(O) 7^ (since V'n(*-') — 0)- Using (j4.3p and the facts > and 
Vb(0) = for the numerical example (12. 6|) . we get 

'^2n-l,2n(^)(^n,V'n) ^ , ^^^^ 2^ 



2V^2(0)c^„ 



Therefore, the bifurcation occurs for 5 > if u}2n~i 2n(^) < (^-S- 
the right of uji) and for 5 < if W2„_i 2n(^) > *-* (^-S- to the left of ujq 

and 104)^ see Table 1. Note that the negative values of 5 correspond to the 
values of 5 below the level 6 = d due to periodicity of the function Gi {6, k) 
in S. The above local existence analysis for even bound states is confirmed 
by the full lines near 5 = in Fig. [5] (b) and near 5 = d = 10 in Fig. [6] (b). 

(ii) Bifurcation of Odd Eigenf unctions 

Similarly to (i) , we study the zeros of G2 (S, k) . We compute the deriva- 
tives of G2{6, k) in 6 and k at {6, k) = (0, 0) 

dsG2(,0,0) = u'_{0) = i^UO), 
9^G2(0,0) = -MO), 

such that the leading-order approximation for the root of G2 {S, k) near 
{5, k) = (0, 0) is given by 

^-^(0)" + ^^"^- 2(C(0))2 - + 0[.), 

where Vn(0) / and D(0) = -^(,(0)^^„(0) (since Vn(0) = 0). Therefore, 
the bifurcation occurs for 5 > if oj2n-i 2n(^) < (^-S- for a; to the right 
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of ws) and for 5 < if W2n~i 2n(^) > (^-S- ^or to to the left of uj2). The 
dashed Unes near 5 = in Fig. [6] (b) and near 5 = d = 10 in Fig. [5] (b) 
confirm this analysis. 

Note that there are curves in Figs. [5](b) and[6](b) which do not bifurcate 
from 6 = and 5 = d = W but still bifurcate from the band edge a; = 
LVn- Bifurcations of these curves cannot be confirmed from the analytical 
theory above, unless the values of Gi^2(^; 0) for < 5 < d are approximated 
numerically. 

5 Conclusion 

We have employed methods of bifurcation theory for the existence problem 
of SGS's supported by the nonlinearity interface and the periodic potential. 
Two bifurcation problems are considered numerically. The first bifurcation 
takes place from the standard gap solitons existing at the zero jump of 
the nonlinearity coefficient. The second bifurcation takes place from the 
bound state consisting of parts of two standard gap solitons glued together 
in a continuously differentiable SGS. Three asymptotic results are described 
in the article. We show that the standard gap solitons can be continued 
generally for small jumps in the nonlinearity coefficient. On the contrary, 
no SGS's for non-zero jump of the nonlinearity coefficient exists in the NLS 
approximation which is valid near the band edges. In addition, we study 
analytically bifurcations of eigenvalues of the Schrodinger operator with non- 
smooth potential from band edges of the Hill's equation. 

One can argue that the SGS's bifurcating from a standard gap soliton 
or a gluing combination of two gap solitons inherits stability properties of 
gap solitons in the neighborhood of the local bifurcation points. Stability of 
standard gap solitons was considered analytically and numerically in [13]. 
The stability properties can change far from the bifurcation points. Detailed 
computatitons of stability of the SGS's will be the subject of the forthcoming 
work. 
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